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SUPERSONIC  WAVE  DRAG  OF  NONLIFTING  DELTA  WINGS  WITH 
LINEARLY  VARYING  THICKNESS  RATIO 
By  Arthur  Henderson,  Jr. 


SUMMARY 


The  supersonic  wave  drag  of  a nonlifting,  symmetrical,  double- 
wedge -profile,  delta  wing  the  thickness  ratio  of  which  varies  linearly 
in  the  spanwise  direction  is  calculated  hy  means  of  linear  theory. 

In  general  it  is  found  that  a delta  wing  with  linearly  varying 
thickness  ratio  can  have  less  wave  drag  than  a constant-thickness-ratio 
delta  wing  of  the  same  plan  form  when  both  wings  have  either  the  same 
projected  frontal  area  or  the  same  internal  volume.  The  thickness 
distributions  for  minimum  drag  and  the  corre sponding  values  of  the  ratio 
of  the  drag  of  a linearly  varying  thickness -ratio  wing  to  a constant- 
thickness -ratio  wing  are  found. 


INTRODUCTION 


In  reference  1,  Puckett  has  found  the  supersonic  wave  drag  of  a 
nonlifting,  symmetrical,  double -wedge -profile,  delta  wing  with  constant 
thickness  ratio.  It  is  shown  that  the  drag  coefficient  for  the  delta 
wing  may  be  reduced  below  the  two-dimensional  value  only  for  the  case 
in  which  both  the  leading  edge  and  the  ridge  line  are  subsonic,  the 
maximum  reduction  for  a given  Mach  number  and  semiapex  angle  being  a 
function  of  the  position  of  the  ridge  line. 

In  the  present  paper  the  additional  effect  on  the  drag  of  varying 
the  thickness  ratio  in  the  spanwise  direction  is  determined  with  the 
assumption  that  the  thickness  ratio  varies  linearly  in  the  spanwise 
direction,  which  means  that,  when  the  wing  is  viewed  from  behind,  the 
line  of  maximum  thickness  has  a parabolic  shape.  The  source -distribution 
method  developed  in  reference  1 is  used  to  represent  the  variable- 
thickness-ratio  wings. 

It  is  shown  that  the  deviation  of  the  maximum-thickness  line  from 
a straight  line  can  be  represented  by  a nondimensional  parameter.  Charts 
are  presented  from  which,  if  this  parameter  is  known,  the  drag  of  a 
variable -thickness -ratio  delta  wing  can  be  found.  Also  shown  is  the  fact 
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that  a delta  wing  with  a linearly  varying  thickness  ratio  can  give  less 
wave  drag  than  a constant -thickness-ratio  delta  wing  when  both  wings  have 
the  sane  projected  frontal  area  or  when  both  have  the  same  internal 
volume.  The  value  of  the  nondimens  ional  maximum-thickness  -distribution 
parameter,  for  which  the  ratio  of  the  drag  of  the  variable -thickness- 
ratio  wing  to  the  drag  of  the  const ant -thickness -ratio  wing  is  a minimum, 
is  calculated  for  both  criteria. 


SYMBOLS 


Af 

a = PeRL 

B =-f- 
creLE 

b = 


projected  frontal  area 


3(y/<=r) 

b 


p = \/yfi  - 1 


Cp  drag  coefficient,  D/qS 

cr  root  chord 


c(y) 

D 

Ap 

eLE 


eRL 


F1>F2.>F3 

G1jG2'G3 

Hi,H2,H3 


>- 


local  chord 
drag 

local  static  pressure  minus  free -stream  static  pressure 
tangent  of  leading-edge  semiapex  angle 
tangent  of  ridge -line  semiapex  angle 

functions  defined  in  appendix  B 

slope  of  airfoil  surface  in  streamwise  direction 
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M 

m 

m 


^T) 

P]_  and  p-^ 


free -stream  Mach  number  • 

dimensional  thickness-distribution  parameter 

nondimens ional  thickness-distribution  parameter, 

mcrgLE  = mcrb 
t (3t 

dummy  variables  for  x and  y,  respectively 

pres sure -distribution  function  and  80q/8x  associated 
with  basic  source  distribution  when  b < 1 


P2  and  p2  pressure-distribution  function  and  associated 

with  superimposed  source  distribution  when  a < 1 

Qp,Q]_  and  q^q-^  pressure-distribution  functions  and  80q/8x  associated 

with  basic  source  distribution  when  b > 1 

Q2,Q2  and  q2,q2  pres sure -distribution  functions  and  802/dx  associated 

with  superimposed  source  distribution  when  a > 1 


<1 


P 

S 

. ^r 
t(y) 


dynamic  pressure. 


free-stream  density 

% 

plan-form  area 
root  thickness 

local  maximum  thickness 
♦ 

thickness  ratio  at  root,  tr/cr 


V 

Vi 

0 


free-stream  velocity 
internal  volume 

perturbation  velocity  potential. 


perturbation  velocity  in  z-direction 
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w 


X-tJt  z 


Subscript 


s: 


Cartesian  coordinates  (right-handed  system) , x in 
direction  of  free -stream  velocity 


i 


E 

1 


right  half  of  wing 
basic  distribution. 


superimposed  distribution 
Primes  denote  constant-thickness-ratio  ving. 


ANALYSIS 


s 

5 

1 

i 


The  linearized  partial  differential  equation  for  the  perturbation 
potential  in  steady  supersonic  flow  is 

i 


p2  _ -JO*  _ _ o 

p ax2  ay2  az2 


(1) 


When  considering  a thin  symmetrical  airfoil  at  zero  angle  of  attack, 
for  which  small  disturbances  are  assumed,  the  boundary  conditions  on 
the  surface  of  the  wing  may  be  satisfied,  to  the  first  order,  in  the 
plane  of  the  wing.  Hence 


> 


(2) 


In  reference  1,  it  is  shown  that  a solution  of  equation  (l).  which  satis- 
fies condition  (2)  is  that  for  a source  distribution 


0(x,y) 


w(6,Tl)d|  dtj 

\A*  - £)2  - P2(y  - n) 


2 > 


(3) 
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where  w(I,tj)  = X(l  ,tj)V  is  proportional  to  the  source  strength  per 
unit  area  and  X is  the  slope  of  the  airfoil  section  measured  in  the 
streamwise  direction.  For  the  wing  shown  in  figure  1,  the  slope  of 
the  surface  ahead  of  the  ridge  line  is  equal  to  the  ratio  of  one -half 
the  local  maximum  thickness  at  the  ridge  line  to  the  part  of  the  local 
chord  ahead  of  the  ridge  linej  the  slope  of  the  rear  surface  is  the 
negative  ratio  of  one-half  of  the  local  maximum  thickness  at  the  ridge 
line  to  the  remainder  of  the  local  chord.  In  the  present  paper  the 
thickness  ratio  is  assumed  to  vary  linearly  with  y,  or 


t(y) 

c(y) 


W 


where  tr/cr  is  the  root  thickness  ratio,  and  m is  a parameter  which 
determines  the  shape  of  the  thickness  curve  in  the  spanwise  direction. 
Therefore,  over  the  forward  part  of  the  wing 


v(x,y) 


t(y)y 

2(1  - r)c(y) 


V / tr 
(1  - r)y2cr 


+ 


(5a) 


and  over  the  rear  part 


w^y) 


= _ t(y)V 
2rc(y) 


(5b) 


where 


r 


fLE 

6RL 


and  represents  the  ridge -line  position  in  percent  chord 


measured  from  the  trailing  edge. 


In  the  following  analysis  30/dx  is  evaluated  from  equation  (3)  hy 
integrating  over  the  right  half  of  the  upper  surface  of  the.  wing.  Inas- 
much , as  the  wing  shape,  and  therefore  the  pressure  distribution,  is  sym- 
metrical with  respect  to  the  x-axis,  50(x,y)/Sx  for  the  entire  upper 
surface  is 
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>vk 


80(x,y) 

8$(x,y) 

+ 

80(x,-y) 

8x 

8x 

R 

8x 

(6) 


A 


The  potential  of  any  point  (x,y)  is  found  by  integrating  over  the  area 
of  the  wing  lying  within  the  forward-drawn  Mach  lines  from  the  point  (x,y) 
in  order  to  obtain  the  effect  on  the  point  of  all  the  sources  that  can 
influence  it. 

The  wing  of  the  present  paper  is  represented  by  a distribution  of 
sources  in  the  plane  z = 0.  This  distribution  is  uniquely  determined 
by  the  boundary  conditions  on  the  wing.  Since  equation  (3)  is  a solu- 
tion to  a linear  partial  differential  equation,  this  distribution  may 
be  represented  by  a superposition  of  source  distributions  in  any  manner 
desired  as  long  as  the  net  result  is  that  prescribed  by  the  boundary 
conditions.  The  drag  on  the  wing  is  then  found  by  computing  the  drag 
caused  by  the  pressure  field  of  each  source  distribution  and  algebrai- 
cally adding  the  results. 

In  the  present  analysis  two  source  distributions  are  used:  one, 

the  basic  distribution  which  extends  over  the  entire  plan  form  and 
follows  the  law  governing  the  source  strengths  representing  the  forward 
slopes  of  the  wing,  and  the  other,  the  superimposed  distribution  which 
lies  in  the  region  bounded  by  the  ridge  lines  and  the  trailing  edge. 

The  sources  in  this  region  have  the  strengths  prescribed  by  the  boundary 
conditions  for  this  region  plus  the  negative  value  of  the  strengths  of 
the  first  distribution.  Thus,  when  superimposed,  the  net  source  distri- 
bution for  the  whole  wing  is  that. prescribed  by  the  boundary  conditions. 

From  figure  2(a),  the  basic  source  distribution  used  to  represent 
a delta  wing  with  subsonic  leading  edge  ( PexiE  ^ 1')  caQ  expressed  as 


R n(l  - r) 


Sy  rx-PCy-q) 


0 Jq/eLE 


^2c, 


+ mT))di  dr) 


\A*  - S)2  - pS(y  - 1)2 


gLE(x+py)  rx+My-T)) 
1+3  e 


n(l  - r) 


LE 


f tr  \ 

2-  + mn  d!  dn 


ti/c^  \/(x  - i)2  - p2(y  - r\)2 


(T) 
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and.  the  superimposed  source  distribution  when  the  ridge  line  is  sub- 
sonic (Psrl  < 1),  as 


gs(x,y)]R  = 


nx“l3(y“T))  / tr  \ 

te + “t p ^ 


”r(1  " ^ !0  Jri/HgL  \/(x  - 5)2  - S2(y  - 1)8 


«r(l  - r) 


i€RL(x+fty) 

1+36rl 


px+p(y-T]) 


U 


2c, 


+ mrijdl  dt] 


n/eRL  vU  - £)2  - 32(y  - b)2 


(8) 

By  the  use  of  equations  (6)  and  (7),  a solution  for  the  basic  distribu- 
tion is  obtained: 


80i 

8x 


2Vte> 


tanh 


-1 


|3jt(l  - r)\/l  - b2  ^ 1 - u2 


2Vmb2x  * VI  - b2  , 

^2  tarLh  + 


P2*(l  - r)(l  - b2)J/  Vl  - u2 


2Vmb->ux 


P2jt(l  - r)(l  - b2) 


3/2 


tanh 


-1  u \Zl~  ~ ft2  + 

b \/l  - u2 


2Vmb2x 


p2jc(l  - r) (l  - b2) 


Vi  - u2 


= Pi(m,r,x,u,b) 


(9) 
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where  u = p— , b = Psle,  and  a = 06rl*  Equation  (9)  is  obtained  in 

X 

appendix  A.  Similarly,  equations  (-6)  and  (8)  produ'ce  the  solution  for 
the  superimposed  distribution: 


ox  p2 


p1(m,r,x,u,a) 


(10) 


The  expression  for  the  potential  of  the  basic  .source  distribution 
used  to  represent  a delta  wing  the  leading  edge  of  which  is  supersonic 
(P6LE  > 1)  is  similar  to  equation  (7);  however,  since  in  this  case  the 
Mach  lines  from  the  apex  lie  on  the  wing,  there  are  two  regions  of  inte- 
gration (fig.  2(b)).  One  is  the  region  interior  to  the  Mach  lines  where 
the  right  and  left  halves  of  the  wing  are  interdependent  and  is 

exactly  equation  (7).  The  integral  is  evaluated  differently,  however, 
since  P^le  > 1.  The  other  is  the  region  that  lies  between  the  Mach 
lines  and  the  leading  edge  where  the  flow  is  two-dimensional.  Here  the 
right  and  left  regions  are  independent  of  each  other,  and  the  expression 
for  the  potential  on  one  side  also  gives  the  potential  for  the  other. 

For  the  region  of  the  wing  exterior  to  the  Mach  lines. 


-P(y-'n) 


1 2c, 


+ mq  d£  dq 


\J(x  - I)2  - P2(y  - q)2 


*(1  - r) 


(x+Py)  rsX+p(y-rj ) 


1+P6LE 


l2c. 


+ mq  d£  dq 


q/e^  \/(x  - I)2  - P2(y  - q)2 


(ID 


which,  with  equation  (7),  yields  for  a supersonic-leading-edge  delta 
wing,  when  0 ^ u ^ 1 
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2V[ — £)b 
V2cr/ 

Ml  - r)  v/b2  - 


, _1  Vb^  - 1 

tan  ± + 


2Vmb2x 

P2n(l  - r ) (b2  - l)3^2 


. -1  \A~  - 1 


rl  - u£ 


2Vmb3ux  . _i  uVb2  - 1 

— tan  x — — ; ; 

- rUb2  _ 1 ^3'2  b \/l  - u2 


32tc(1  - r ) (b2  - I)' 


2Vmb2x 


(32n(l  - r ) (b2  - l) 


a - u£ 


= q^a^r^u^) 


and,  when  1 = u ^ b 


¥fe> 


3 

Vmb  ux 


Vmb2x 


3(1  - r) \/b2  - 1 32(1  - r)(b2  - l)  ' 32U  - r)(b2  - 1) 


= q-,  (m,r,x,u,b) 


The  corresponding  functions  for  the  superimposed  source  distribu 
tion  when  the  ridge  line  is  supersonic  are  for  0 ^ u ^ 1 


= qo  = 


- z q-,  (m,r,x,u,a) 
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and,  for  1 ^ u ^ a. 


&r = q2 


- ^ qx(m,r,x,u,a) 


(13b) 


The  parameter  m in  the  preceding  equations  has  the  dimensions  of 

(length)^".  It  will  he  useful  to  examine  the  curves  of  thickness  plotted 
against  spanwise  distance  that  can  he  produced  hy  a variation  of  m 
(fig.  3)  with  the  purpose  of  forming  a new  .dimensionless  parameter  which 
includes  m.  From  equation  ( 4) , a nondimens ional  form  of  the  thickness 


curve  is 


iM . C(y/i  + fffi 

tr  lcr  tr 


and  since,  from  figure  1,  c(y)  = cr  - 


iM  . 1 + (?*  . _i_V  . 2e_ 

tr  \ creLE/  treLE 


= 1 + 


2mcreLE 


^mcrgLE  ( j \2 


creLE 


T \creLE 


= 1 + 


2mcrb  Wy/er)  2mc  h pfy/c 


(15a) 


where 


The  nondimens ional  thickness -distribution  parameter  is,  therefore,  taken 
to  be  mcrh/(3T  which  is  hereinafter  denoted  by  m.  Equation  (15a)  may 
be  rewritten  as  follows: 
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11 


= 1 + (2m  - l)B  - 2mB2  (15b) 

lr 


where 

■ B _ g(y/cr) 

b 

Equation  (15b)  is  a quadratic  in  B and  represents  a parabola.  If  the 
type  of  curve  for  equation  (15b)  is  known,  one  absolute  limit  on  5 can 
be  established.  This  limit  can  be  determined  since,  at  the  wing  tip  the 
slope  of  the  thickness  curve  cannot  be  less  than  zero.  For  this  case, 

m = -Tj.  When  the  slope  of  the  thickness  curve  is  zero  at  the  root, 

m = however,  this  value  is  not  necessarily  an  upper  limit.  Curves  of 

thickness  distribution  for  different  values  of  m are  shown  in  figure  3- 


PRESSURE  DISTRIBUTION 


To  the  first  order,  the  pressure  on  the  surface  of  the  wing  is 
given  by 


(16) 


Introducing  the  parameter  5 into  equations  (9),  (10),  (12a),  (12b), 
(13a),  and  (13b)  and  using  equation  (l6)  gives  the  pressure  equations  in 
the  following  forms: 


For -the  basic  subsonic-leading -edge  delta  wing 


&E. 

9 


2Tb 


tanh 


0*(1  - r)\Jl  - b2 

Ml  - r)(l  -'b2)U  _ b2 


-1  Q 

VT  - u2 


tanh 


-1  Vi  - 12 

\/l  - u^ 


_hL_  t^-1  ™/LzE  . 

1 - b2  bVl  - U2  / 


= P(m,r,^,u,b^  = Px 


(0  £ u < b)  (17) 
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For  the  superimposed  subsonic -leading-edge  delta  wing 

^ p(m,r ;^;u,aj  = P2  (0  ^ u < a)  (18) 

For  the  basic  supersonic-leading -edge  delta  wing 


= Q^m,r,^pu,b^  = Qx  (l  ^ u < b)  (19^) 


For  the  superimposed  supersonic-leading-edge  delta  wing 


If  = = "?  §2 


(0  ^ u 5 1) 


^ = -r  Q/m.r .-r-,u ,a^  = -r  Qo  (1  S u £ a) 


(20b) 
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The  pressure  field  associated  with  the  basic  source  distribution 
is  shown  in  figure  k for  both  the  subsonic-  and  supersonic -leading -edge 

cases  for  particular  wings.  For  these  wings  m = which  prescribes 

that  the  maximum-thickness  line,  which  is  at  the  trailing  edge,  have 
maximum  concavity  (see  fig.  3)*  The  effect  of  having  m < 0 is  that, 
as  x/cr  increases,  the  pressures  are  progressively  decreased  in  a 
spanwise  direction  below  the  values  for  a conical  pressure  field.  If 
m > 0,  the  pressures  would  be  progressively  increased  in  a spanwise 
direction  above  the  conical  pressures  as  x/cr  increases. 


DRAG 


For  the  drag  of  the  delta  wing  herein  considered,  the  following 
three  distinct  cases  must  be  treated: . 

Case  I - Supersonic  leading  edge  and  supersonic  ridge  line  (b  > 1,  a > 1) 
Case  II  - Subsonic  leading  edge  and  supersonic  ridge  line  (b  < 1,  a > 1) 
Case  III  - Subsonic  leading  edge  and  subsonic  ridge  line  (b  < 1,  a < 1) 


In  computing  the  drag  by  the  method  of  superposition,  it  must  be 
kept  in  mind  that  the  basic  source  distribution  creates  a pressure  field 
the  influence  of  which  is  felt  over  the  whole  surface  of  the  wing; 
whereas  the  pressure  field  from  the  superimposed  source  distribution  is 
felt  only  over  that  part  of  the  wing  which  lies  between  the  ridge  lines 
and  the  trailing  edge,  except  for  case  III  where  the  pressure  field 
"spills  over"  the  ridge  lines  and  exerts  its  influence  on  the  forward 
slopes  up  to  the  Mach  lines  from  the  ridge -line  apex. 

Since  the  wings  are  symmetrical  about  both  the  x,y-  and  x,z-planes, 
the  drag  is  computed  for  the  right  half  of  the  upper  surface  of  the  wing 
and  the  result  is  then  multiplied  by  four. 


Case  I.-  The  drag  for  case  I (fig.  5(a))  is  represented  by 


° = r^x  ds 
<1  Js  * 


(21) 


which,  as  a result  of  the  superposition  of  source  distributions,  becomes 
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h.  ~ / Q1X1  63  + / %>2  dS  “ 

J Ax  ^B1+B2 


QlX2  dS  - ± 


Q2X<2  dS 


The  limits  of  integration  on  the  second  and  fifth  integrals  make 
the  operation  rather  cumbersome.  This  difficulty  may  be  obviated  by 
rewriting  equation  (22)  in  the  following  form,  which  includes  the  sub- 
1 - r 

stitution  X2  X-j_: 


_D_  _ 1 
tap  r 


= | f as  - i-fT  r as 

v/A^"riBp+B2 


as  + i f %Xj. 

^b- 


dS  - 


— ■-  F ds  + -----  r 

r j A2+B3  r2  J b2+B3 


Q2A.^  dS 


For  this  case  b > 1,  a > 1. 

The  element  of  area  is  dS  = dx  dy,  but  the  variables  of  integration 
are  x and  u where  u = 3^.  Therefore,  if  y = u£  and  the  Jacobian 
of  x and  y is  taken  with  respect  to-  x and  u,  the  element  of  area 

X 

becomes  dS  = -^  dx  du.  Performing  the  operations  indicated  in  equa- 
tion (23),  reducing  the  result  to  coefficient  form,  and  making  the  sub- 
stitution a = — yields 
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cDP 


= F]_(r,b)  + 5F2(r,b)  + 5fTo(r,b)  (b  > 1 > r)  (2k) 


Equations  for  F-j_,  F2,  and  F^  are  given  in  appendix  B and  are 

plotted  in  figure  6. 

Case  II.-  The  drag  for  case  II  (fig.  5(B))  is  given  by 


1 j“  = T pl^l  as  + £ P^  dS  - ^ f QgXg  dS  - i f dS 


,b1+b2 


Bn 


Br 


- P^!  dS  - 


1 - ’r 


/A+Bq+B2 


P^  dS  + --  ~ -r  / Q/2k1  dS  + 

r2  JBl 


dS 

r JB0 


(25) 


In  evaluating  equation  (25)  it  must  be  remembered  that  b <1,  a > 1. 
The  final  value  may  be  expressed  in  coefficient  form  by 


cDP 


= G1(r,b)  + mG2(r,b)  + m2G^(r^b) 


(1  > b > r)  (26) 


Equations  for  G^,  G2,  and  G3  are  given  in  appendix  B and  plotted 

in  figure  7 with  the  H-functions  which  are  presented  subsequently. 

Case  III.-  The  drag  for  case  III  (fig.  5(c))  is  expressed  as 


_D 


= r P1x1  dS  + f p^g  dS  - - f P2X2  dS  - - f P 
Ja1+A2  ^B  T r JAo 


dS 


- I 

r 


f hh.  is  - 

Ja„  • - 


JAj+ Ag 


L 


P X-l  dS  + ^ f PgXi 

A1+A2+B  . r2  J3 


- fA  p2^i  as 


dS  - 


(27) 


1 6 


MCA  TN  2858 


Here  b < 1,  a < 1^  and  the  result-  is 

= Hr(r,b)  + mH2(r,b)  + 52H3(r,b)  (l  > r > b)  (28) 

T2 


Equations  for  H^,  H2,  and  H^  are  given  in  appendix  B and  plotted 

in  figure  7- 


Whether  a variable-thickness-ratio  delta  wing  is  better  from  the 
wave-drag  standpoint  than  a constant-thickness-ratio  delta  wing  depends 
upon  the  criterion  chosen.  Equations  (2k) , (26),  and  (28)  give  the  drag 
for  a const ant -thickness -ratio  delta  wing  when  m = Oj  the  results  are 
identical  to  those  of  reference  1.  Inasmuch  as  F2  > F3,  G2  > G3, 


and  H2  > H3  always , it  is  obvious  that,  if  the  variable-  and  constant- 
thickness -ratio  wings  are  compared  on  the  basis  of  identical  plan  form, 
r,b,  and  thickness  ratio  at  the  root,  the  concave  thickness-distribution 

curve  /-i  ^ m < 0,  see  fig.  3^  will  give  lower  drag.  This  criterion. 


however,  is  a poor  one  because  the  internal  volume  and  the  projected 
frontal  area  are  reduced.  Since  the  same  projected  frontal  area  appears 
to  provide  a fair  basis  for  comparison  and  internal  volume  is  important 
in  practical  design  considerations,  these  two  criteria  are  developed  for 
the  purpose  of  comparing  the  drag  of  a constant-  with  a variable -thickness 
ratio  delta  wing. 


PROJECTED -FRONTAL -AREA  CRITERION 


If  a variable -thickness -ratio  delta  wing  and  a constant -thickness- 
ratio  delta  wing  have  identical  plan  form,  r,  b,  and  projected  frontal 
area,  then 


Af 


2 


1 creLE 


t(y)dy 


and  from  equation  (15b) 


t(y)  = tr 


2m  - 1 

y 

creLE 


2m  2\ 
cr2SLE2  7 ) 
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Therefore , 


(29) 


where  Af  is  the  projected  frontal  area,  S is  the  plan-form  area, 
and  t is  the  thickness  ratio  at  the  root. 


For  a constant-thickness -ratio  delta  wing 


Af'  = S*t' 


(30) 


From  equations  (29)  and  (30),  the  stipulation  that  Af  = Af ' and 
S = S'  gives  the  relationship  between  r'  and  r as 


r 1 


From  equation  (24) 


CD13  _2 

— = Fx  + mF2  + m F3 

"r 

for  the  variable-thickness-ratio  delta  wing,  and 


Cp’3 


- F1 


(3D 


(32) 


for  the  constant -thickness-ratio  delta  wing.  Forming  the  ratio  of  equa- 
tion (24)  to  equation  (32)  and  substituting  equation  (31)  gives: 


F-j_  + mF2  + m^F3 

(1  + f s)\ 


(33) 
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The  ratio  Cp/Cp'  can  he  minimized  with  respect  to  m,  and  the  value 
of  m,  which  gives  minimum  drag  for  the  conditions  stipulated,  is 


4^  - 3F2 

6f3  - 2F2 


(34) 


Equation  (34)  applies  to  the  G-  and  H-functions  as  well  as  to  the 
F-functions  where  the  corresponding  G-  and  H-functions  are  substituted 
for  the  F-functions. 

The  nondimens ional  thickness -distribution  parameter  m for  minimum 
drag,  based  on  frontal -area  considerations,  is  presented  in  figure  8. 
Figure  9 presents  Cp/Cp 1 as  calculated  from  equation  (33)  using  the 
values  of  m obtained  from  figure  8.  When  b > 1 > r,  it  can  be  seen 
from  figure  8 that  the  wing  for  minimum  drag  is  essentially  that  with 
constant  thickness  ratio  (m  % 0),  except  when  b = 1.001,  and  even  for 
this  case  figure  9 shows  a maximum  drag  reduction  of  less  than  2 percent 

When  1 > b > r,  figure  8 indicates  that  the  wing  for  minimum  drag 
has  a concave  thickness  distribution  (m  < 0).  The  corresponding  curves 
of  figure  9 indicate  a maximum  drag  reduction  of  only  about  4 percent 
for  this  case.  The  main  conclusion  to  be  drawn  then  is  that,  with  the 
projected  frontal  area  held  constant,  the  thickness  can  be  shifted  from 
the  tip  to  the  root  with  no  penalty  in  wave  drag. 

When  1 > r > b,  the  optimum  thickness  distribution  is  sometimes 
concave  (m  < 0,  for  which  the  same  conclusion  can  be  drawn  as  for 
1 >b  > r)  and  sometimes  convex  (m  >0).  Some  of  the  indicated  shapes 

for  m >0  axe  highly  unrealistic  ^ that- is,  m » which  indicates 

that,  within  the  realm  of  physically  practical  shapes,  there  is  no 
minimum  m.  This  condition,  however,  does  not  mean  that  a drag  reduc- 
tion is  not  possible.  Figure  10  presents  the  curve  of  Cp/Cp'  plotted 
against  m for  a case  when  m,  as  determined  from  equation  (34),  gives 
an  unrealistic  shape.  It  can  be  seen  that,  although  m for  minimum 
drag  as  calculated  from  equation  (34)  is  3-68,  values  of  m from  l/2 
to  1 give  a drag  reduction  of  20  to  25  percent.  In  figure  9 the  dashed 
lines  show  values  of  Cp/Cp * when  m represents  a wing  that  is  unreal- 
istic, where  an  unrealistic  wing  has  arbitrarily  been  chosen  to  be  one 
for  which  m > 1.  (See  fig.  3.) 
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INTERNAL -VOLUME  CRITERION 


In  order  to  compare  the  variable-  and  constant-thickness-ratio 
delta  wings  on  the  basis  of  internal  volume , both  wings  are  specified 
to  have  identical  plan  form,  r,  b,  and  internal  volume.  The  internal 
volume  of  the  wing  of  figure  1 is 


Vi  = 


z dx  dy 


(35) 


For  the  right  half  of  the  wing,  from  the  leading  edge  to  the  ridge  line 


(36) 


and  from  the  ridge  line  to  the  trailing  edge 


z 


t(y) 

c 


lfl(Cr  ~X) 


(37) 


Using  equations  (35) } (36),  and  (37)  gives,  for  a variable-thickness- 
ratio  wing. 


V. 

1 


(38) 


and  for  a const ant -thickness -ratio  wing 


t’c. 


3 


(39) 


When  equations  (38)  and  (39)  are  equated  and  the  conditions  that 
Vi  = Vi’,  cr  = cr',  and  S = S'  are  fulfilled. 
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t ' = t (1  + ^ m 


(to) 


CDp  Cd'3 

Forming  the  ratio  of  to  , as  before  for  the  present  criterion. 


i2 


gives 


n > 

LD 


Fjl  + 5F2  + 


1 -A2 

2 / 


1 + £ m)  Fx 


(41) 


which,  when  minimized  with  respect  to  m,  gives 


5 = 2 


F1  - F2 
4f3  - F2 


(42) 


for  the  value  of  m for  minimum  drag  under  the  conditions  stipulated. 

Values  of  5 for  minimum  drag  based  on  internal -volume  considera- 
tions are  presented  in  figure  11.  Figure  12  presents  Cp/Cp ’ as  calcu- 
lated from  equation  (4l)  by  the  use  of  the  values  of  m from  figure  11. 
When  b > 1 > r,  figure  11(b)  shows  that  the  value  of  5 for  minimum 
drag  is  approximately  constant  at  about  m = -0.45;  that  is,  the  thickness 
distribution  for  minimum  drag  is  very  close  to  that  of  maximum  concavity. 
The  corresponding  curve  of  figure  12  gives  a maximum  drag  reduction  of 
between  8 and  10  percent. 

When  1 > b > r,  figure  11(a)  indicates  a value  of  5 which  is 
fictitious  in  that  m which  would  prescribe  a wing  where  part  of 

the  upper  surface  is  below  the  bottom  surface  and  conversely.  Therefore, 
there  is  no  minimum  drag  for  any  real  wings  in  this  range.  Figure  13  is 
a plot  of  Cp/Cp 1 against  m for  a particular  configuration  in  this 

range.  Although  the  optimum  5 as  determined  from  equation  (42)  is 
fictitious,  it  can  be  seen  that  any  value  of.  5 between  5 = and 

5=0  will  give  a drag  reduction,  the  reduction  in  drag  at  m = -~ 

being  about  20  percent.  In  figure  12,  the  values  of  Cp/Cp ' which  were 
obtained  by, using  fictitious  values  of  5 are  represented  by  short 
dashed  lines. 
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When  1 > r > b,  figure  11  shows  that  the  m for  minimum  drag 
varies  from  negative  to  positive  with  increasing  r for  each  value  of  b. 

Included  in  this  range  are  both  fictitious  ^m  < -7^  and  unrealistic' 

^m  » 7^  values  of  m calculated  from  equation  (42).  The  corresponding 
values  of  Cp/Cp 1 are  shown  in  figure  12  where  the  short  dashed  lines 
were  obtained  by  the  use  of  fictitious  values  of  3 ^ and  the 

long  dashed  lines  indicate  the  use  of  unrealistic  values  of  m (m  >1). 

It  should  be  emphasized,  however,  that  substantial  dra‘g  decreases  can  be 
realized  in  the  regions  of  the  dashed  curves  for  practical  wings,  as 
can  be  seen  from  figures  10  and  1-3,  which  indicate  the  trend  of  all  curves 
of  Cp/Cp'  plotted  against  m when  the  m for  minimum  Cp/Cp'  is 
either  positive  or  negative. 


ILLUSTRATIVE  EXAMPLE 


If,  for  a given  constant-thickness-ratio  delta  wing,  it  is  desired 
to  increase  the  thickness  of  the  root  and,  for  structural  reasons,  the 
thickness  at  the  tip  is  not  to  be  decreased,  two  methods  are  suggested. 
Either  a variable-thickness-ratio  wing  with  prescribed  root  and  tip  thick- 
ness ratios  or  a new  constant-thickness-ratio  wing  with  prescribed  root 
thickness  ratio  could  be  formed. 

For  the  former  case,  if  t ' is  the  thickness  ratio  of  the  original 
const ant -thickness -ratio  wing,  the  slope  of  the  maximum  thickness  line 
with  respect  to  y is  -T'/2eLE-  If  the  variable -thickness-ratio  wing 
has  this  slope  at  the  tip,  equation  (15)  gives 


&L 


8 

¥l  2 Jv_r  g 2eTF 

y^creLE  ^ 


T (1  + 2m)  = - T 


2e 


LE 


or 


t'  = t(1  + -2m) 


(43) 


If  the  root  thickness  of  the  original  wing  is  increased  by  n times. 
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t = riT1  (MO 

The  two  equations  (43)  and  (44)  yield  the  necessary  thickness  distribu- 
tion 


m 


1 - n 
2n 


(45) 


As  an  example , consider  the  wings  for  which  b = 0.8  and  r = 0.5.  If 
the  root  thickness  of  the  original  constant-thickness-ratio  wing  is 

increased  by  l/2,  n = ^ and  5 = When  the  ratio 

<=■  6 


Cjq  Gq_  + mGg  + m^G^ 

CD*  (1  + 25)^ 


is  formed,  the  drag  of  the  variable -thickness -ratio  wing  is  found  to  be 
1.77  times  that  of  the  original  constant-thickness-ratio  wing. 

If  the  increased  root  thickness  is  attained  by  forming  a new 
constant-thickness-ratio  wing,  the  drag  is  2.25  times  that  of  the 
original  wing.  Here,  of  course,  both  the  projected  frontal  area  and 
the  internal  volume  axe  larger  than  those  for  the  variable -thickness- 
ratio  wing.  The  drag  of  a constant-thickness-ratio  wing  which  has  the  . 
same  projected  frontal  area  as  the  previously  mentioned  variable-thickness- 
ratio  wing  is  I.78  times  the  drag  of  the  original  wing;  whereas,  the 
drag  of  a constant-thickness -ratio  wing  which  has  the  same  internal 
volume  as  the  previously  mentioned  variable-thickness-ratio  wing  is  I.89 
times  that  of  the  original  wing. 


CONCLUDING  REMARKS 


The  variable -thickness -ratio  delta  wing  has  been  compared  with  the 
constant-thickness-ratio  delta  wing  under  the  conditions  that  they  have 
identical  projected  frontal  area  and  identical  internal  volume.  For 
these  conditions  the  optimum  value  of  the  thickness-distribution  parameter 
has  been  determined.  The  wing  shape  given  by  this  value,  although  not 
always  practical  or  real,  nonetheless  gives  an  indication  of  the  types 
of  maximum-thickness-line  distribution;  that  is,  whether  it  should  be 
convex,  straight,  or  concave.  From  the  analysis  the  following  conclusions 
are  drawn: 
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1.  On  the  basis  of  fixed  frontal  area, 

(a)  When  the  leading  edge  and  ridge  line  are  both  supersonic, 
the  constant -thickness -ratio  delta  wing  has  essentially  optimum 
wave-drag  characteristics  for  all  combinations  of  leading -edge 
semiapex  angle  and  position  of  the  ridge  line. 

(b)  When  the  leading  edge  is  subsonic  and  the  ridge  line  is 
supersonic,  the  concave  maximum -thickness  distribution  has  the 
optimum  wave-drag  characteristics  which,  compared  with  a constant- 
thickness -ratio  wing,  gives  greater  thickness  at  the  root  and  less 
thickness  at  the  tip.  It  should  be  noted,  however,  that  the  drag 
reduction  which  can  be  realized  in  this  range  is  small. 

(c)  When  the  leading  edge  and  ridge  line  are  both  subsonic, 
as  the  ridge  line  moves  forward  from  the  sonic  condition  to  the 
leading  edge,  the  maximum-thickness  distribution  for  optimum  wave 
drag  goes  from  concave  to  straight  to  convex.  For  the  latter  case, 
the  drag  decrease  indicated  is  large  and  absurd  values  are  approached 
as  the  ridge  line  approaches  the  leading  edge. 

2.  On  the  basis  of  internal  volume  considerations, 

(a)  When  the  leading  edge  and  ridge  line  are  both  supersonic, 
the  concave  maximum-thickness  distribution  has  optimum  wave-drag 
characteristics,  the  drag  reduction  being  about  8 to  10  percent  for 
all  combinations  of  leading-edge  semiapex  angle  and  position  of  the 
ridge  line. 

(b)  When  the  leading  edge  is  subsonic  and  the  ridge  line  is 
supersonic,  the  concave  thickness  distribution  has  the  best  wave- 
drag  characteristics,  with  drag  reductions  of  as  much  as  20  percent 
being  predicted. 

(c)  When  the  leading  edge  and  ridge  line  are  both  subsonic, 
as  the  ridge  line  moves  forward  from  the  sonic  condition  to  the 
leading  edge,  the  maximum-thickness  line  goes  from  concave  to 
straight  to  convex;  whereas  the  corresponding  drag  reductions  range 
from  about  20  percent  to  zero  to  absurd  values. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  October  1,  1952. 


2k 


NACA  TN  2858 


Let 


Then 


EVALUATION  OF 
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Differentiating  (0l)^  with  respect  to  x gives 


By  the  use  of  equations  (l6o)  and  (169)  in  reference  2,  the  following 
equation  is  obtained: 
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(a)  Case  I:  b > 1,  a > 1. 


\ 

(b)  Case  II:  b < 1,  a > 1.  (c)  Case  III:  b < 1,  a < ] 

Figure  5*-  Regions  of  integration  for  computing  the  drag. 


(c)  F3(r,b). 

Figure  6.~  Concluded. 


(a)  Gi(r,b)  and  Hi(r,b).  G-functions 
are  on  the  right  of  the  cusp. 

Figure  7 . - Variation  of  G-  and  H-functions  with  1 - r ’ for  different 

values  of  b. 


(b)  'and  B^r^b).  G-functions 

are  on  the  right  of  the  cusp. 

Figure  7«-  Continued. 
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(c)  63(2*, t>)  and  H3(r,b).  G-functions 

are  on  the  right  of  the  cusp. 


Figure  7.-  Concluded. 


(a)  m for  G-  and  H-f unctions. 


ft) 


5 for  F-functions, 


Figure  8.-  Variation  of  m for  minimum 
frontal-area  considerations  with  1 - 


Cp/Cp'  based  on  pro jected- 
r for  different  values  of  b. 
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